Computing the distance between quantum channels: 
Usefulness of the Fano representation 
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The diamond norm measures the distance between two quantum channels. From an operational 
vewpoint, this norm measures how well we can distinguish between two channels by applying them 
to input states of arbitrarily large dimensions. In this paper, we show that the diamond norm can 
be conveniently and in a physically transparent way computed by means of a Monte-Carlo algorithm 
based on the Fano representation of quantum states and quantum operations. The effectiveness of 
this algorithm is illustrated for several singie-qubit quantum channels. 

PACS numbers: 03.65.Yz, 03.67.-a 



I. INTRODUCTION 

Quantum information processes in a noisy environment 
are conveniently described in terms of quantum chan- 
nels, that is, linear, trace preserving, completely positive 
maps on the set of quantum states [l], Q • The problem 
of discriminating quantum channels is of great interest. 
For instance, knowing the correct noise model might pro- 
vide useful information to devise efficient error-correcting 
strategies, both in the fields of quantum communication 
and quantum computation. 

It is therefore natural to consider distances between 
quantum channels, that is to say, we would like to quan- 
tify how similarly two channels £i and £2 act on quantum 
states, or in other words to determine if there are input 
states p on which the two channels produce output states 
p'l = £i{p) and p'2 = £2{p) that are distinguishable. The 
trace norm of p'l — p'2 represents how well p'l and p'2 can be 
distinguished by a measurement ^ : the more orthogonal 
two quantum states are, the easier it is to discriminate 
them. The trace distance of two quantum channels is 
then obtained after maximizing the trace norm of p'^ — p'2 
over the input state p. However, the trace norm is not 
a good measure of the distance between quantum chan- 
nels. Indeed, in general the presence in the input state of 
entanglement with an ancillary system can help discrim- 
inating quantum channels \^]^ ■ This fact is captured by 
the diamond norm the trace distance between the 

overall output states (including the ancillary system) is 
optimized over all possible input states, including those 
entangled. 

The computation of the diamond norm is not known 
to be straightforward and only a limited number of al- 
gorithms have been proposed [Tll - [l3| . based on compli- 
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cated semidefinite programming or convex optimization. 
On the other hand, analytical solutions are limited to 
special classes of channels 0, II]. In this paper, we pro- 
pose a simple and easily parallelizable Monte-Carlo al- 
gorithm based on the Fano representation of quantum 
states and quantum operations. We show that our al- 
gorithm provides reliable results for the case, most sig- 
nificant for present-day implementations of quantum in- 
formation processing, of single-qubit quantum channels. 
Furthermore, in the Fano representation quantum opera- 
tions are described by affine maps whose matrix elements 
have precise physical meaning: They are directly related 
to the evolution of the expectation values of the system's 
polarization measurements [ll, [3 [TH] . 

The paper is organized as follows. After reviewing in 
Sec.|TT]basic definitions of the distance between quantum 
channels, we discuss in Sec. IIIII two numerical Monte- 
Carlo strategies for computing the diamond norm. The 
first one is based on the Kitaev's characterization of the 
diamond norm. The second one is based on the Fano rep- 
resentation of quantum states and quantum operations. 
The two methods are compared in Sec. lIVI for a few phys- 
ically significant single-qubit quantum channels. Finally, 
our conclusions are drawn in Sec. |Vl 



II. THE DIAMOND NORM 

A. Basic definitions 

We consider the following problem: given two quantum 
channels £1 and £2, and a single channel use, chosen uni- 
formly at random from {£1, £^2}, we wish to maximize the 
probability of correctly identifying the quantum channel. 
It seems natural to reformulate the optimization problem 
into the problem of finding an input state (density ma- 
trix) p in the Hilbert space Ji such that the error prob- 
ability in the discrimination of the output states £i{p) 
and £2{p) is minimal. In this case, the minimal error 
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probability reads 



/ _ 1 ll'S'l — ^2||l 

~ 2 4 ' 

||f 111 = max||£'(p)||i, 

p 



(1) 



where ||X||i = TryXtX denotes the trace norm. 

The superoperator trace distance ||£i — £2||i is, how- 
ever, not a good definition of the distance between two 
quantum operations. The point is that in general it is 
possible to exploit quantum entanglement to increase the 
distinguishability of two quantum channels. In this case, 
an ancillary Hilbert space /C is introduced, the input state 
^ is a density matrix in /C (8) 'H, and the quantum opera- 
tions are trivially extended to fC. That is to say, the out- 
put states to discriminate are {1^. ® f i)C and {Ik: ® ^2)^, 
where I/c is the identity map acting on /C. The minimal 
error probability reads 



1 11^1-^2 



Pe 



||£-|U^max||(ZK®£)e||i, 



(2) 



where \\£\\o denotes the diamond norm oi £. It is clear 
from definition ^ that 



\\£\U = \\Tk®£\\i > mil 



(3) 



and therefore pe "£ p'ej so that it can be convenient to 
use an ancillary system to better discriminate two quan- 
tum operations after a single channel use. The two quan- 
tum channels £1 and £2 become perfectly distinguishable 
{pe — 0) when their diamond distance ||£i — £2\\o = 2. 

It turns out that the diamond norm does not depend 
on /C, provided dim(A^) > dim(H) Q. Due to the con- 
vexity of the trace norm, it can be shown that the maxi- 
mum in both Eqs. ([T]) and ^ is achieved for pure input 
states 161. 



Note that the space H is traced out in the definition of 
^A, rather than the space TZ. Finally, it turns out 
that [i [l3 



ll^llo =i^max(*A,*s), 



(6) 



where Fmax(^'A, ^'s) is the maximum output fidelity of 
and ^_B, defined as 

i^max(*A,*B) =maxi^[*A(/5l),*s(p2)], (7) 

P1.P2 

where pi,p2 are density matrices in and the fidelity 
F is defined as 



*b) = Tr^J ^Y'^-B^y. (8) 

Note that are not density matrices: the condi- 

tions Tr(*A) = 1, Tr(\l/B) = 1 are not satisfied. 



III. COMPUTING THE DIAMOND NORM 

We numerically compute the distance (induced by the 
diamond norm) between two quantum channels £i and £2 
using two Monte-Carlo algorithms. The first one is based 
on the direct computation of | \£i — ^21 |o, with the output 
states {Xic £i)S, and {Xjc £"2)^ in Eq. ^ computed 
from ^ taking advantage of the Fano representation of 
quantum states and quantum operations. The second 
Monte-Carlo algorithm uses the Kitaev's representation 
of the diamond norm to compute the maximum output 
fidelity Fmax of Eq. ([7]). In the following, we will refer 
to the two above Monte-Carlo algorithms as F-algorithm 
and K-algorithm, respectively. For the sake of simplicity 
we will confine ourselves to one-qubit quantum channels, 
even though the two algorithms can be easily formulated 
for two- or many-qubit channels. 



A. The F-algorithm 



B. Kitaev's characterization of the diamond norm 

Kitaev provided a different equivalent characterization 
of the diamond norm, see, e.g., [I,[l3|- Any superoperator 
(not necessarily completely positive) £ : L{H) — )■ LCH), 
with L{'H) space of linear operators mapping % to itself, 
can be expressed as 

£{X) - Tvn{AXB^), (4) 

where X € L{T-L), A and B linear operators from % to 
TZ ® T-L, with TZ auxiliary Hilbert space and dim (7?.) < 
[dim('H)]. It is then possible to define completely positive 
superoperators Ai'^ b ■ L{Ti) L{TZ): 

^A{X)^TTn{AXA^), *b(^) -Tr„(BXBt). (5) 



In this section we describe the F-algorithm, which we 
will use to directly compute the diamond norm with 
the maximum taken over a large number of randomly 
chosen input states ^. A convexity argument shows 
that it is sufficient to optimize over pure input states 
^ = |5')(^'| 16]. For one-qubit channels, it is enough to 
add a single ancillary qubit when computing the diamond 
norm [9] . Therefore, we can write 

I*) = CoolOO) -I- CoilOl) -f CiollO) + Ciilll), (9) 

with 

Coo = cos 6*1 cos 612, 

Coi =cos0isin02e"^S 

Cio = sin6llCos6'3e*'^^ ^ ' 

Cii = sin 01 sin 6*36*, 
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where the angles 6i G [O, -l] and the phase (f)i e [0, 2tt]. 
Hence, the maximization in the diamond norm is over the 
6 real parameters ^1,^2, ^3, and 4'i:4'2,4'3- Of course, 
the number of parameters can be reduced for specific 
channels when there are symmetries. 

Let £1 and £2 denote the two single-qubit superoper- 
ators we would like to distinguish. The output states 
^[ = {X/c (8) £1)^ and ^2 = i^ic ^ £2)^ can be conveniently 
computed using the Fano representation. Any t wo-q ubit 
state can be written in the Fano form as follows [17Hl9| : 



Ra/BCTa <8) (7/3, 



(11) 



-x,y,z,I 



where a^;, fXy, and are the Pauh matrices, aj — 1, and 

= Tr[(a„®a^)^]. (12) 

Note that the normalization condition Tr(^) = 1 implies 
Rij = 1. Moreover, the coefficients RajS are real due to 
the hermicity of ^. Eqs. (fTTj) and (fT2|) allow us to go from 
the standard representation for the density matrix (in the 
{|0) = |00), |1) = |01), |2> EE 1 10), |3) EE 111)} basis) to the 
Fano representation, and vice versa. It is convenient to 
write the coefRcients R^p as a column vector, 

-f^ [Rxx: Rxy-j Rxz-j Rxl-j Ryx: Ryy : 
Ryz J Ryl 7 Rzx 7 Rzy j Rzz j Rzl j Rjx ; Riy ; Rlz ; Rll] 

(13) 

Then the quantum operations andljc'8'f2 map, in 

the Fano representation, R into R'^ = Ai\ H. and R'a ~ 



A^2 R-) respectively, where M.Y'' and M.2''' are affine 
transformation matrices. Such matrices have a simple 
block structure: 



(2) 



<(2) 



(14) 



with X^^'' and 4x4 affine transformation matri- 

ces corresponding to the quantum operations Xjc and 
£i (of course, X^^^ is the identity matrix). Matrices 

Aif'^ directly determines the transformation, induced by 
£i, of the single-qubit Bloch-sphere coordinates {x,y,x). 
Given the Fano representation for a single qubit, p = 
\Y.a'^a(Ta, then r = [x^y^z^iy and y[ = A^f^r. We 
point out that, while one could compute 4' from the 
Kraus representation of the superoperator I;^ 'S^ £ii the 
advantage of the Fano representation is that the matrix 
elements of Mi are directly related to the transforma- 
tion of the expectation values of the system's polarization 
measurements [H. [l^. [T5| . 

Finally, we compute the trace distance between ^[ and 

as 



where Ai, A4 are the eigenvalues of ^[ — £,2- 



(15) 



B. The K-algorithm 

We consider a special and unnormalized state in the 
extended Hilbert space 1C®'H: 
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|j/c)|jw), 



(16) 



where dim(/C) = dim(H) and {Ijk;)}, {|jw)} are or- 
thonormal bases for /C, T-L. The state \a) is, up to a nor- 
malization factor, a maximally entangled state in IC^Ii. 
Wc define an operator a on /C (g) "H : 



cj = {XK,®£){\a){a\), 



(17) 



where £ — £1 ~ £2 '\s the difference of two quantum oper- 
ations but is not a quantum operation itself. That is, £ is 
linear but not trace preserving and completely positive. 

Using the singular value decomposition of ti, we can 
write 



M 



a = ^|«W)(z; 



(18) 



where M is the rank of a and ju^'-*), Iw*-*-*) are vectors in 
JC ®T-L. The operator a completely specifies £ and can 
be exploited to express £ as 



M 



(19) 



The operators A^^'^ and B*^'^ can be derived by gener- 
alizing the construction of the operator-sum representa- 
tion for quantum operations (see, for instance, Sec. 8.2.4 
in Ref. Q]). Given a generic state IV') = V'jli'w) in T~L, 

we define a corresponding state in /C: |-0) — V'JIjk;)- 
Next, we define 



For instance, in the single-qubit case 



,112 , U3 , U4 



and we obtain 



U2 u\' 



(20) 



(21) 



Similarly, we define 



(22) 



Finally, it can be checked that with the above defined op- 
erators A''^\ B^'^'^ we can express £ by means of Eq. (|19p . 

We can now give explicit expressions for a and b 
in Eq. ©: 



(23) 
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(24) 



with 1 < i, j < M. Therefore, and are M x M 
matrices. In the single-qubit case, M < A and to compute 
the diamond norm through Eq. ([7|) we need to calculate 
eigenvalues and eigenvectors of matrices of size M. A 
simpler but less efficient decomposition off, is 
provided in Appendix |^ 

For single-qubit channels the optimization ([7]) is over 
6 real parameters, 3 for pi and 3 for p2 (for instance, 
the Bloch-sphere coordinates of pi and p2)- As dis- 
cussed in Sec. IIIIBl the same number of parameters are 
needed in the F-method. However, the F-method has 
computational advantages in that only the eigenvalues of 
the 4x4 matrix — ^2 are required, while in the K- 
method we need to evaluate both eigenvalues and eigen- 
vectors of matrices in general of the same size {M — 4). 
Moreover, the singular-values decomposition of matrix 
a is needed. Besides computational advantages, the F- 
method is physically more transparent, as it is based on 
affine maps, whose matrix elements have physical mean- 
ing, being directly related to the transformation of the 
expectation values of the system's polarization measure- 
ments 



IV. EXAMPLES FOR SINGLE-QUBIT 
QUANTUM CHANNELS 

In this section, we illustrate the working of the F- and 
K-methods for the case, most significant for present-day 
implementations, of single-qubit quantum channels. 



A. Pauli channels 

We start by considering the case of Pauli channels, 

cx—x,y,z,I OL 

(25) 

for which the diamond norm can be evaluated analyti- 
cally 17|: 



(26) 



this value of the diamond norm being achieved for maxi- 
mally entangled input states. The Pauli-channel case will 
serve as a testing ground for the F- and K-algorithms and 
help us develop a physical and geometrical intuition. 

Let us focus on a couple of significant examples. We 
first consider the bit-flip and the phase-flip channels: 



^AP) = — —(^xP^^x-, 

1 + C2 ^ — Cl 

^2[P) = — ^o-^pc^z, 



(27) 
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FIG. 1: Schematic drawing of the trace (and diamond) dis- 
tance distance between the bit-flip and the phase-flip chan- 
nels. 



with < ci, C2 < 1. We then readily obtain from Eq. ((26 
that 



\\£-\\o = 11^1 - '^^2||o = max{l - ci, 1 - C2}. 



(28) 



For this example, the diamond norm coincides with the 
trace norm, ||f||o = II^Hi, and therefore has a simple 
geometrical interpretation: For a single qubit the trace 
distance between two single-qubit states is equal to the 
Euclidean distance between them on the Bloch ball 
Superoperators E\ and £2 map the Bloch sphere into an 
ellipsoid with x (for the bit-flip channel) and z (for the 
phase-ffip channel) as symmetry axis. If we call (x, y, z), 
(x']^, yj., -^i)' and (3^21^21-^2)1 the initial Bloch-sphere co- 
ordinates and the new coordinates after application of 
quantum operations £\ and £2, respectively, we obtain 



ClZ, 



(29) 



Z/2 = C2?;, 



The geometrical meaning of the trace norm for the 
present example is clear from Fig. [T] the length of the 
line segment pp'2 — p'\p'i is the trace (and the diamond) 
distance ||£i— 52II1 = ||i^^i~^2||o (note that in this figure 

1 - C2 > 1 - Cl). 

As a further example, we discuss a special instance of 
the channels considered in Ref. Q: 



^ , , 1 1 1 

^Ap) =2^"^ -^f^xPOx + ^O-yPO-y, 

'^^2(p) = o^pa^. 



(30) 



In this case, 



(31) 
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FIG. 2: Schematic drawing of the trace distance distance be- 
tween the two Pauli channels of Eq. (|30p . 



The trace norm ||£^||i = |, as shown in Fig. [21 where 
||£||i is given by the length of the segment p[p'2- On 
the other hand, the two channels £i and £2 are perfectly 
distinguishable, as we readily obtain from Eq. (|26l) that 
||£||o = 2, this value being achieved by means of maxi- 
mally entangled input states. Therefore, in this example 
||f||o > ||£^||i, that is, entangled input states improve the 
distinguishability of the two channels. 

Channels ((30)) are very convenient to illustrate the con- 
vergence properties of the F- and K-algorithms. Let 
us first consider the F-algorithm. If the maximum of 
IICi-'f2lli' with 4' = {lic<S)£t)£. {i = 1,2) is taken over A^^ 
random initial conditions, then the obtained value differs 
from the diamond norm by an amount S{Nr) which 
must converge to zero when Nr — > 00. To get conver- 
gence to \ \£\\o for channels (|30p it is enough to optimize 
over real initial conditions. Numerical results, shown in 
Fig. [3] for a few runs, are consistent with S{Nr) ^ l/Nr- 
A rough argument can be used to explain the l/Nr de- 
pendence. Assuming a smooth quadratic dependence of 
the distance D{£^) = — C2II1 the parameters for 
optimization [the angles 9i in (fTO|) . with 9i S [Oj f ] j 
6*2,^2 S [0,27r)] for ^ around the value ^0 optimizing D, 
then |D(^) — I?(^o)| ^ S when the distance between ^ 
and ^0 is of the order of v^- The number Nr of ran- 
domly distributed initial conditions typically requested 
to get a point satisfying ||C — Colli < VS is of the order 
of {1/VS)"'', where Up is the number of parameters for 
optimization. Therefore, S{Nr) ^ (l/Nr)'^/^. In the ex- 
ample of Fig. H np = 3 leading to S{Nr) - (l/Nr)^/^. 
On the other hand, numerical data exhibit a l/Nr depen- 
dence. This fact has a simple explanation: the maximum 
is achieved for Bell states, which are invariant under ro- 
tations. Hence, the maximum distance is obtained not 
on a single point but on a curve and the number Nr of 



00 




FIG. 3: Error S in evaluating the diamond distance between 
channels (|30p by means of the F-algorithm after A'^,. randomly 
chosen initial conditions. Three different runs are shown, to- 
gether with the l/Nr (fuU line) and {1/Nrf^^ (dashed line) 
dependences. 

initial conditions requested to get a point within a trace 
distance smaller than ^/S from this curve is of the order 
of {l/Vsy'''-\ thus leading to 6{Nr) - l/Nr. 

With regard to the K-algorithm, we have observed in 
the example of Eq. ([50)) the same 1 /Nr convergence to the 
expected asymptotic value \\£\\o = 2. However, the cost 
per initial condition in the K-algorithm is much larger 
than in the F-algorithm, in agreement with the general 
discussion of Sec. lIIIBl Furthermore, the physical mean- 
ing of the F- method is much more transparent. For in- 
stance, for the Pauli channels ([25)) the affine transforma- 
tion matrix A^f ^ such that R; = {l'-^^ <^ mI^^)R. has a 
simple diagonal structure: 

= diag[(c,)„ (cy),, (c,),, 1], (32) 

where 

{Cx)i = 1 - 2[{qy)i + (q^),], 

{cy), = 1 - 2[iq,), + (q,),], (33) 

iCz)i = 1 - '2[iqx)t + iqy)i]- 

Matrix M^/'^ simply accounts for the transformation, in- 
duced by £i, of the polarization measurements for the 
system qubit: a[ — {ca)ia, with a — x,y, z Bloch-sphere 
coordinates. 

We have checked the computational advantages of the 
F-method also for all the other examples discussed in this 
paper. For this reason in what follows we shall focus on 
this method only. 

B. Nonunital channels 

In this section, we consider nonunital channels, that is, 
channels that do not preserve the identity. Therefore, in 



71/2 



FIG. 4: Trace norm II^Hi — \\£i — falli, where £i and £2 
are displacement channels along the +x- and +2-axis of the 
Bloch sphere. 



contrast to the Pauli channels considered in Sec. IIV A[ 
single-qubit nonunital channels displace the center of the 
Bloch sphere. Such channel are physically relevant in the 
description, for instance, of energy dissipation in open 
quantum systems, the simplest case being the amplitude 
damping channel jli, i2l . 

As a first illustrative example, we compute the dis- 
tance between two channels £1 and £2 corresponding to 
displacements of the Bloch sphere along the +x- and +z- 
direction ij] . In the Fano representation the affine trans- 
formation matrices A^|^'' corresponding to maps have 
a simple structure: 



MY 



M 



(1) 











^x 



















Cx 














1 










■ 



















CI 


SI 











1 



(34) 



(35) 



where we have used the shorthand notation Cx = cos Ox 
and Cz = cos Oz- The two channels depend parametri- 
cally on Bx^dz G [O, |] . The limiting cases 9x — ^ and 
0z = Q correspond to f 1 — I and £2 = I, respectively. 
Moreover, for 9x — ^ quantum operation £1 maps the 
Bloch ball onto the single point [x'l — l-.y'i — 0, z[ — 0); 



for Oz 



the mapping operated by £2 is onto the north 



pole of the Bloch sphere, {x'2 — 0,y2 — 0, 23 — 1). 

The numerically computed trace distance \\£i — £2\\i 
is shown in Fig. 21 as a function of the parameters 6x 
and 9z- We gathered numerical evidence that for such 
channels the diamond distance equals the trace distance. 

It is interesting to examine the analytical solutions for 
two limiting cases: (i) 6x = (the case 6*^ = is anal- 



FIG. 5: Bloch-sphere visualization of the trace norm ||f ||i in 
the limiting cases 6x ~ (left) and 9x = 6z = (right). 



ogous) and (ii) 0^ = = f . In the first case, £1 maps 
the Bloch sphere into an ellipsoid with z as symmetry 
axis while £2 = I- The trace norm is given by the 
length of the line segment = p'\p'2 shown in Fig. [5] 
(left). 



|£||i = 2sin2 



(36) 



It is interesting to remark that, in contrast to the Pauli 
channels, the optimal input state is not a maximally en- 
tangled input state. In the limiting case (i) the channels 
are strictly better discriminated by means of an appro- 
priate separable input state, i.e., the south pole of the 
Bloch sphere (see the left plot of Fig. [5])) rather than by 
Bell states. Indeed, given the maximally entangled input 
state 



1 



!*)(*!, |*> = _(|oo) + |ii)), 



(37) 



we obtain 



Cz - 1 


52 
a-1 -si 



(38) 



The trace distance between i^x ^2 is then computed 
by means of Eq. 03 



Sl-^St^^(\-Czf 



(39) 



As shown in Fig.[i H^J - C2II1 < \\£\\\ for any Cz > 0. 

In case (ii), for any initial state p, the Bloch coordi- 
nates of p'l and p'2 are given by (1,0, 0) and (0, 0, 1). The 
trace norm is given by the distance between these two 
points, \ \£\\i = [see Fig. [5] (right)]. In this case, given 
an input Bell state or any other two-qubit input state ^, 
the trace distance — C2II1 — — \\£\\i- Thus, there 
is no advantage in using an ancillary system. 
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0.1 0.2 0.3 0.4 0.5 



FIG. 6: Comparison between \\£i — £2\\i (upper curve) and 
(lower curve) for = f . 



As a further example, we compute the distance be- 
tween the depolarizing channel £i and the nonunital 
channel £2 corresponding to the displacement along the 
-|-z-axis of the Bloch sphere. The depolarizing channel 
belongs to the class of Pauli channels with qj = 1—p 
and = Qy = Qz — ^- This channel contracts the Bloch 
sphere by a factor (l — with < p < |. 

Fig.[7]shows the numerically computed \\£\\i and ||£||o 
as well as their difference \\£\\o — \ \£\\i- The use of an 
ancillary qubit improves the distinguishability of the two 
channels. However, we remark again that maximally en- 
tangled input states can be detrimental. For instance, 
in the limiting case p = 0, 9^ = ^ we obtain from 

Eqs. dSS]) and jSSl) ||£||i = 2 > - ^2111 = A 
clear advantage is instead seen in another limiting case, 
p = J, 9z = 0- The fully depolarizing channel maps 
each point of the Bloch ball onto its center, so that the 
trace distance is given by the radius of the Bloch sphere, 
I — f 2 1 1 1 = 1 • On the other hand, by means of Eq. ([26|) 
we obtain H^i — f2|U = §■ 



V. CONCLUSIONS 

We have shown that the distance between two quan- 
tum channels can be conveniently computed by means of 
a Monte-Carlo algorithm based on the Fano representa- 
tion. The effectiveness of this algorithm is illustrated in 
the case, most relevant for present-day implementations 
of quantum information processing, of single-qubit quan- 
tum channels. A main computational advantage of this 
algorithms is that it is easily parallelizable. Furthermore, 
being based on the Fano representation, is enlights the 
physical meaning of the involved quantum channels: the 
matrix elements of the affine map representing a quan- 
tum channel directly account for the evolution of the ex- 
pectation values of the system's polarization measure- 
ments. More generally, we believe that the Fano rep- 
resentation provides a computationally convenient and 
physically trasparent representation of quantum noise. 




FIG. 7: From top to bottom: trace norm ||f||i — \ \£i —£-2\\i, 
diamond norm ||f ||o, and their difference ||f ||o — ||f ||i; £1 is 
the depolarizing channel, £2 the displacement along the +z- 
axis of the Bloch sphere. 



Appendix A: Alternative decomposition of £ 

Given a superoperator £ = £1 — £2, with £1 , £2 quan- 
tum operations, we start from the Kraus representa- 
tion of f 1 and £2- 

Ah M2 

£,{X) = Y,EkXEl f2(A) = ^FfcAF,^ (Al) 

k=l k=l 
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and define new operators: 

_g(2fc-l) ^ ^(2fc)^ j^(2k) ^ ^(2fc-l)^ 

_' ' (A2) 

where fc = 1, M = max(Mi, M2). Note that, if Mi > 
M2, we set Gfe = for k = M2 + 1, Mi; vice- versa, if 
Ml < A/2, Fk = for fc = Afi + 1, ...,M2. It is easy to 
see that 



quired, but less efficient. Indeed, the maximum num- 
ber of terms in (jA3l) is twice that of decomposition (|19l) . 
This implies that, if 5*^, VPs are expressed in terms of 
operators rather than A^'), S^'), to evaluate 

F{^a,^b) we need to compute eigenvalues and eigen- 
vectors of matrices of size 2M. In the single-qubit case, 
typically 2M = 8. 



2M 



£{X) = ^i^XB 



wt 



(A3) 
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